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We show the details of the pathway for dissociative ionization process of ground electronic state 
of aligned due to first ionization of D2 in short (~100 fs) and intense (4.0 x 10^"^ W cm~^) 
480 nm laser pulses. The initial vibrational state of D J comes from the vertical transformation of 
the ground state of D2. The initial wavepacket in the ground electronic state of is outgoing 
through dissociation-ionization channel accompanied by a strong coupling between IsGg and 2pau 
electronic states. We show explicitly that the transition from the coupling states IsGg and 2pau 
to the ionization state is not a direct transition but takes place through other intermediate states 
with some dissociation energy that results in the internuclear distribution of the ionization to move 
considerably to larger internuclear distances. 



PACS numbers: 33.80.Rv, 33.80.Gj, 42.50.Hz 
I. INTRODUCTION 

When molecules are exposed to an intense laser pulse 
(even for the simplest molecules, i.e. H2 and H2"), a 
wealth of fascinating phenomena may be observed. A 
few examples are: above-threshold or tunneling ioniza- 
tion, charge resonance enhanced ionization, dissociative- 
ionization, above threshold dissociation, bond softening 
and hardening, high order harmonic generation, etc. [l]. 
An interesting subject in studying of H2 (D2) in the in- 
tense laser pulse is double ionization [2 . The double 
ionization of H2 (D2) exposed to an intense laser field 
can proceed through either non-sequential or sequential 
mechanisms, depending on the intense laser pulse param- 
eters. Sequential double ionization (SDI) is the dominant 
mechanism at higher intensities and also for circularly 
polarized fields. 

Another interesting new phenomenon has been ob- 
served in a recent experimental study of D2 exposed to an 
intense laser field that revealed a new high-energy band 
for the kinetic energy release (KER) spectrum at short 
wavelength in the SDI [3] . In the experiments of Litvinuk 
et. al., the linearly polarized laser pulses of different peak 
intensities with duration of ~100 fs were focused into a 
well-collimated supersonic jet of target D2 molecules in- 
side the uniform-electric- field ion imaging spectrometer. 
They found a new high-energy band for the KER spec- 
trum in an 480 nm laser pulse with peak intensity 4 x 10^^ 
W cm~^. The emergence of this new band depend on the 
wavelength, i.e. the energy of this band is decreasing with 
increasing the wavelength until this band merges with the 
enhanced ionization band for 800 nm and longer wave- 
length [3 . The position of this high-energy band seems 
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FIG. 1: We suppose that at a definite time, Ton, H2 (D2) 
is ionized vertically from its ground state to the ground elec- 
tronic state of (DJ)- Ionization of (DJ) practically 
occurs through the DIG. The nuclear component in this chan- 
nel (DIG) posses both dissociation energy (DE) and Goulomb 
explosion energy (GEE). The green arrows are related to the 
results of the present simulation. 



to be independent of the intensity for the various wave- 
lengths. 

To explain this phenomenon, Litvinyuk and et. al. 
proposed some possible pathways to emerge of this high- 
band energy (Fig.jl]) as follow. After the neutral molecule 
is singly ionized, the vibrational wavepacket propagates 
on the field-dressed electronic potential energy surfaces 
of D2" , staying mostly in the lowest two electronic states 



2 



(field- free dg and states) . All pathways involve three- 
photon resonance of the ground state to the first ex- 
cited electronic state of D2", Lq. In one pathway, Li 
m Fig. [ij the ion's internuclear distance increases in a 
dissociative state until reaching the critical values {Rc)^ 
where the second ionization takes place with high proba- 
bility and the Coulomb explosion occurs. The final KER 
of the nuclear fragments will include the kinetic energy 
gained during dissociation and also kinetic energy of the 
Coulomb explosion {l/Rc). In another possible pathway, 
L2 in Fig.[l] the second ionization may take place directly 
at shorter internuclear distances for the shorter wave- 
lengths. In this case the three-photon excitation to the 
(Ju state is only a virtual one, with ionization yield being 
strongly enhanced by the presence of the intermediate 
resonance. The fragments would gain their full kinetic 
energy on the Coulomb potential, and that energy will 
be very close to that of the first discussed pathway, Li, 
since the Coulomb and potentials are nearly parallel 
to each other. The present experimental data are not suf- 
ficient to determine which of these two pathways may be 
the real pathway or near to reality [3 . Another possible 
pathway, a net-two-photon pathway (three-photon exci- 
tation followed by subsequent de-excitation at the one- 
photon resonance), would result in lower kinetic energy 
fragments that is inconsistent with the experimental re- 
sults [3 , 8 . 

To explore the details of the Litvinyuk's experiment 
and the real pathway, it is necessary to set up a sim- 
ulation that allows us to declare the details of the 
dissociation- ionization pathway of D2". The dynamics 
of II2 (D2) and H2" (D2") exposed to an intense laser 
field are very complicated because it involves electron 
and nuclear dynamics simultaneously with different time 
scales that results in a complex dissociation-ionization 
process. In intense laser field, electron dynamics occur 
in attosecond time scales and nuclear dynamics, i.e. vi- 
bration and rotation, takes place in femtosecond and pi- 
cosecond time scales respectively. It is sometimes possi- 
ble based on Born Oppenheimer approximation (BOA) 
to investigate the nuclei and electron dynamics, sepa- 
rately. This approach is extensively used to investigate 
electronic dynamics of molecules in intense laser fields 
pT| [2m I2T] . But this assumption fails when both nuclei 
and electron dynamics are important and considerable. 
This situation happens with intense ~ 10- 100 femtosec- 
ond laser pulses, when the duration of the laser pulse 
is comparable to the vibrational period of molecules [3]- 
[8]. In these conditions, a complex simulation based on 
the solving of the time dependent Schrodinger equation 
(TDSE) beyond the BOA would be required for the com- 
plete description of molecular dynamics. For a molecule 
with two or more electrons, this task is very far from 
present available computer ability, even without consid- 
eration of the nuclear dynamics. The rigorous approach 
is feasible only for linear molecules aligned with the elec- 
tric laser field with only one electron. In this condition, 
it is possible to consider the full dimensional electron dy- 



namics and also the vibrational nuclear dynamics in the 
solution of the TDSE for H+ (D+) 0110]. 

In this work, we have performed a simulation for D2" 
beyond the BOA and also without the soft-core approxi- 
mation [Hi [12], i.e., by the rigorous solution of the TDSE 
for the full dimensional electron dynamics of the aligned 
D2" with the electric laser field. This paper is organized 
as follows: the details of the numerical implementation 
and the simulation setup are presented in Sec. II. Results 
and discussions appear in Sec. III. Section IV contains 
observations and concluding remarks. Throughout this 
article we use the atomic units unless stated otherwise. 



II. NUMERICAL SOLUTION OF THE TDSE 

Time dependent Schrodinger equation in the cylindri- 
cal polar coordinate system for II2" (D2") located in the 
laser field parallel to the internuclear axis in atomic unit 
(au) can be expressed as 

. di,{z,p,R,t) ^ ^^^^ ^ ^^^^^^ 

where the Hamiltonian for this system, H{z^ p^R^t)^ is 
given by [IMS] 
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with £^0 being the laser peak amplitude, u = 27ru the an- 
gular frequency, and finally f{t) the laser pulse envelope 
which is set as 



f{t) exp 



-2\n{2){t-tony 



(4) 



where Tp is the full width at half maximum (FWHM) 
duration of the Gaussian shape of the laser pulse, with 
Tp =25 cycles in this work. For the time discretization 
of the TDSE, a propagator derived from split-operator 
methods has been used. This propagator is unitary and 
is obtained by combining the classical split operator and 
the Crank-Nicholson method [16] . The details of our cal- 
culations are described in our previous reports [3 [151 HZj • 
In this simulation, the time step is set to St = 0.02 and 
5tR = 0.2 for the electron and nuclei time propagation 
respectively and rp=25 cycles (~40 fs). The differen- 
tial operators in Eq. ^ are discretized by the eleven- 
point difference formulae which have tenth-order accu- 
racies [17 . To solve the above TDSE numerically, we 
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adopted a general non-linear coordinate transformation 
for both electronic and nuclear coordinates. For the spa- 
tial discretization, we have constructed a finite difference 
scheme with a non-uniform (adaptive) grid for z and p 
electronic coordinates which are finest near the nuclei and 
coarsest at the border regions of the simulation box. A 
finite difference scheme with an adaptive grid is used also 
for R coordinates that is finest for small R and become 
a coarse grid for large R. Using a fine grid for electronic 
coordinate and p) near the nuclei and for small val- 
ues of internuclear distance coordinate {K) improves the 
treatment of the electron dynamics near the nuclei (the 
Coulomb singularities) and bound states of the nuclear 
dynamics, while the use of a coarse grid near the borders 
improves the speed of calculations. The grid points for 

p, and R coordinates in our simulation are 560, 155, 
and 240 respectively. The finest grid sizes values in this 
adaptive grid schemes are 0.13, 0.2, and 0.025 respec- 
tively for p, and R coordinates. The grids extend up 

^max ^ ^max ^ 25, and i?^^^ = 16. 

The simulation setup is shown schematically in Fig. [l] 
The laser pulse in this simulation suddenly turns on 
at time Ton- We are assuming Ton is the time on 
which a H2" (D2") is suddenly created according a Frank- 
Condon transition from the neutral H2 (D2) to H2" 
(D2") (Jg state (Fig. hj). In an intense laser field, some 
parts of D2" wavepackage become unbound and outgo- 
ing through different channels. A part of this unbound 
wavepacket becomes outgoing as + D through dis- 
sociation channel (DC) and another part becomes out- 
going through dissociation-ionization channel (DIC) as 
-\- (Fig. [1]). The nuclear components in the 
dissociation-ionization channel possess both dissociation 
and Coulomb explosion energies [9\. 



III. RESULTS AND DISCUSSION 

Figure |2] shows the simulation KER spectra and com- 
pares with the experimental result. In this figure, the 
dotted red line is related to the experimental result [3 for 
the 480 nm linearly polarized laser pulse with 1=4.0 x 10^^ 
W cm~^ intensity and ~100 fs pulse duration. The black 
solid, grey solid, and blue dot-dash lines show simula- 
tion results that are respectively related to 4.0 x 10^^, 
3.5 X 10^^, and 3.0 x 10^^ W cm~^ intensities and a 
480 nm wavelength. The simulation results show ex- 
plicitly two bands, a high-energy and a lo-energy band 
which also have appeared in the experimental results 
|3 . Some inconsistencies can be seen in Fig. [2] between 
the experimental and simulation KER spectra related to 
1=4.0 X 10^^ W cm~^ intensity. These inconsistencies are 
due to some experimental uncertainties in the shape and 
duration of the laser pulse and the unavoidable simplifica- 
tions and approximations applied in the simulation setup. 
However, Fig. [2] shows many similarities and correspon- 
dence between experimental and simulation results. The 
simulation results, in correspondence to the experimen- 



tal results [3], show that the decreasing of the intensity of 
the laser pulse causes a decrease in the signal of the high- 
energy band with respect to the low-energy band (Fig.[2|. 
To derive the details of the KER structure, we represent 
a time-dependent construction of this signal in Fig. [3] As 
can be seen in this figure, at first the high-energy band 
appears between ^4.8 to ^8.1 cycles of the laser pulse 
(note that the meaning of t = in our simulation, as 
shows in Fig.[l] is 2.5 cycles earlier than of the peak of the 
laser pulse). After this duration, it seems that the path- 
way of the ionization with high-energy fragment is closed 
and another pathway of the ionization with low-energy 
fragment is opened. Therefore, this figure proposes that 
there are two distinct pathways for the ionization with 
different energies and time domains. 

Figure |4] represents the space-time population of the 
ground state (Iscr^), first exited state (2pcr^), other non- 
ionization states (ONIS), and total non- ionization states 
(TNIS). The difference between the population of TNIS 
and the population of Isag + 2pau equals to the space- 
time population of the ONIS. This Figure shows explic- 
itly the charge resonance among Isa^, 2p(Ju^ and ONIS. 
The resonance between Iscig and 2pau is very strong and 
a weaker coupling of these two states with the ONIS are 
visible in the ONIS's figure. These figures show explicitly 
the space-time dissociation-ionization pathways of D2". 
Note that, the reduction of the population of the TNIS 
in the dissociation pathways is due to the ionization. Fig- 
ure |4] shows that the population of 2pau becomes consid- 
erable over ^2-5 internuclear distances and ^3-8 cycles 
and for ONIS over ^3-7 internuclear distances and ~5-ll 
cycles. D2" in the dissociative pathways is considerably 
ionized and will has little population before it dissociates, 
as has shown in the next figures. 

To see the details of the dissociation ionization path- 
ways, it is important to derive separately the space- 
dependent population and time-dependent population of 
different states. Fig. [5] represents the time-dependent 
population of the ground state (Isa^, GS), first exited 
state (2p(Tn, ES), the ground and first exited states to- 
gether (GES), other non- ionized states (ONIS), total 
non-ionized states (TNIS) or norm, and finally the time- 
dependent population of the ionized states (IS). As can 
be seen in this figure, the population of the ONIS in- 
creases considerably from about 4 cycles. The charge res- 
onance between different states is observable. A strong 
charge resonance between Gs and Es appears after 2 cy- 
cles results in the population of 2pau becoming more than 
that of IscTg state occasionally. In this figure a weaker 
charge resonance is also observed between the GES and 
ONIS. 

Figure [6] shows the internuclear populations of the ion- 
ization of D2" related to the 4.0 x 10^^ W cm~^ intensity 
in Fig. [2] The internuclear populations of the ioniza- 
tion of D2" in Fi g. [6| ha ve been derived by the virtual 
detector method |15[ I18j . Indeed, this figure represents 
the R-dependent ionization rate, with an arbitrary unit, 
of beyond BOA [9 . The BOA approach has often 
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FIG. 2: Comparison of the experimental and simulation KER spectra. The red dotted line is related to the experimental 
results [3] for a 480 nm laser pulse with peak intensity I==4.0 x lO^'^ W cm~^ and ~100±10 fs pulse duration. The black solid, 
grey solid, and blue dot-dash lines show simulation results for 480 nm, Tp ^40 fs laser pulses with peak intensities 4.0 x 10^^, 
3.5 X 10^^ ,and 3.0 x 10^^ W cm-^ respectively. 
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FIG. 3: Time-dependent construction of the KER Spectrum that is related to the peak intensity 4.0 x 10^^ W cm~^ shown in 
Fig. |2] The high-energy band of the KER Spectrum appears at first between ~4.8 to ~8 cycles of the laser pulse. After this 
time, the low-energy band of the KER Spectrum appears. 



been used to derive the R-dependent ionization rate of 
D2" in most previous studies [2QH22] . In this figure, we 
also show the internuclear populations for the different 
time intervals shown in Fig. [3j This figure shows explic- 
itly that the high-energy band of the KER spectrum in 
Fig. |3j that is constructed between ~ 4.8 — 8 cycles, has 
a strong ionization rate about ~ 4 — 5 internuclear dis- 
tances and the low-energy band of the KER spectrum 



(Fig. |3|, that is constructed between ~ 8.1 — 12.5 cycles, 
has a relatively strong ionization rate about ^6 — 7 inter- 
nuclear distances. After 12.5 cycles, the intensity of the 
laser field becomes considerably low (Fig. [T]) and a signif- 
icant population of D2" was becomes ionized as shown in 
Fig. |5j Therefore it is impossible to appear an enhanced 
ionization signal for larger internuclear distances. 

In Coulomb explosion imaging [23 , we want to re- 
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FIG. 4: Space-time population of the ground state (Iscr^), first exited state (2pau), other no n- ionized states (ONIS), and total 
non- ionized state (TNIS). In this figure, R shows the internuclear distance of DJ. 
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FIG. 5: Time-dependent population of different states of related to the 4.0 x 10^^ W cm ^ intensity in Fig. [2] 
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late the KER spectra to the internuclear population of 
the wavepacket of molecules and so derive the time- 
dependent structure of the molecules. Figures |3]and |6] 
can relate the structure of the KER spectrum and the 
structure of the internuclear population of the ionization. 
If we show the relation of the time-dependent internuclear 
population of the non-ionized D2" wavepacket and the 
time-dependent internuclear populations of the ionized 
D2", it clearly displays the relation of the time-dependent 
structure of molecules and the KER spectrum. In Fig. [7| 
we present the time-averaged internuclear population of 
different states of D2" for the same time intervals that 
is indicated in Fig. [6) The different states in Fig. [7| are 
chosen as follow: the TNIS (solid grey lines), the GES 
(dashed red lines), and the ONIS (dot-dash black lines). 
In this figure, for better comparison, we also show the 
internuclear populations of the ionization in Fig. |6] as IS 
(ionized states) by solid blue lines. Among the differ- 
ent states of D2" in Fig. [rj it is obvious that the ONIS 
is very similar to the IS. For the better comparison of 
the IS and ONIS spectra, in each figure we have scaled 
the overall population of the IS to the overall population 
of the ONIS and presented as scaled ionized state (SIS, 
dot-dot-dash blue lines) in Fig. [7j This figure shows that 
the structure of the SIS is very similar to the structure 
of ONIS and therefore we can conclude that the origin of 
the KER spectrum is from the ONIS and the KER spec- 
trum can relatively map the time-dependent structure of 
the ONIS. 

Figure [3] to [7| constitute a complete set that enable us 
to trace the detail events occurring during dissociative 
ionization of D2" which can be expressed as follows. Dur- 
ing the to ~ 5 cycles. Fig. |5] shows a strong coupling 
between GS and ES and the GS is outgoing through a DC 
(Fig.[4|. In this duration. Fig. [s] shows that the majority 
of the population of the TNIS is in the GES and the pop- 
ulations of the other states are negligible. Near the end 



of this time interval, the population of the ONIS becomes 
considerable. Among different populations in Fig. [7]^b), 
the IS (SIS) and ONIS are very similar but there is a lit- 
tle displacement in the population of the IS with respect 
to the TNIS. This figure also shows that the structure of 
the IS and the GES are similar but the IS shifts to right, 
i.e., to larger internuclear distances with respect to the 
GES. This displacement shows a considerable DE for the 
IS. Fig.[7|^b) also shows that the nuclear fragments of the 
IS have considerable GEE because the IS distribution has 
small internuclear distance distribution. 

During the 5-8 cycles of the laser pulse, the high- 
energy band of the KER spectrum is constructed (Fig.[3|. 
Fig. [rj^b-c) show the internuclear populations of the dif- 
ferent states for this band. Fig. [7]^b-c) also show that 
among different states, the ONIS has the closest struc- 
ture to IS (SIS). Therefore, this confirms that the IS does 
not originate directly from the GES but it arises from the 
ONIS. This point is represented schematically in Fig. [l] 
Generally, all sections in Fig. [7| approve that the struc- 
ture of the IS is very similar to the ONIS, for low or 
high-energy band of KER spectrum. 

Comparing of Fig. [7|^a-c) reveals that during ioniza- 
tion, different states are outgoing through the DG be- 
cause their internuclear populations shift to larger inter- 
nuclear values. In Fig.[7| we can see that from 4.8 to 12.5 
cycles, the magnitude of the displacement of the IS with 
respect to the GES and ONIS is reduced, i.e. the DE of 
the IS is reduced gradually. We can also see that the peak 
of the IS spectrum is placed on the right of the GES spec- 
trum in Fig.jrj^b) opposite to that in Fig.[7|^f). Therefore, 
Fig. [7] shows explicitly that the transition from the GES 
to the IS is not a vertical transition and the population 
transfer from the GES to the IS occurs through the in- 
termediate states (ONIS). Correspondence between the 
IS (SIS) and ONIS spectra in Fig. [T^d-f) are very good 
and the displacement between the GES and the IS (or 
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FIG. 7: IS (solid blue lines) shows outgoing internuclear population from the dissociation ionization channel of for the 
specified time interval on each figure. Other lines as follows: the time-average internuclear population of the TNIS (solid grey 
lines), GES (dashed red lines) and ONIS (dot-dash black lines). To better comparison of the IS and ONIS values, the IS overall 
population are scaled to the ONIS overall population and showed by the dot-dot-dash blue lines (SIS). 



ONIS and the IS) becomes less with respect to Fig. ^h- 
c). Therefore, for these space- time intervals that the 
low-energy band of the KER spectrum is constructed, 
i.e. Fig. [7|^d-f), during the transition of the population 
from the GES to the IS, the ionized fragments obtain the 
smaller DE and their energy come mainly from the CEE. 



IV. CONCLUSION 

In summary, the pathway for dissociative ionization 
process of the aligned ground electronic state of D2" ex- 
posed to a short (~100 fs) and intense (4.0 x 10^^ W 
cm~^) 480 nm laser pulse that its initial vibrational state 
comes from vertical transformation of the ground state 
of D2 can be expressed as follows. During the ~5 cy- 



cles of the laser pulse, the ground state of D2 is strongly 
coupled with 2pau (Fig. [l] and Fig. |5|. Within 3 cycles, 
from ^5 to ~8 cycles of the laser pulse, the high-energy 
band of the KER spectra is constructed (Fig.[3|. During 
this time period, a considerable population of the GES 
transfers to the ONIS and then to the IS (Fig. |5) Fig. [6) 
and Fig. [7|^b-c)). The ionized fragments related to this 
high-energy band is originated mainly from the ~ 4 — 5 
internuclear distances (Fig.jij Fig.[6| and Fig.[7|^b-c)) and 
their internuclear distributions (SIS) are very similar to 
the internuclear distributions of the ONIS, but with some 
small displacements to larger internuclear distances and 
also there is some deformation with respect to the ONIS 
(Fig. [7|^b-c)). Therefore, the structure of the IS in Fig. [t] 
does not map the GES but it is similar to the struc- 
ture of the ONIS. In another words, the transition from 
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the GES to the IS is not direct but this occurs through 
the intermediate states (ONIS) with a considerable DE 
that results in the internuclear distribution of the ioniza- 
tion moves considerably to larger internuclear distances. 
Therefore, an R-dependent displacement between the IS, 
ONIS, and the GES structures are observed. 

After ~8 cycles of the laser pulse, the laser field is still 
strong enough and also the population of the TNIS is 
stih significant. The TNIS (i.e. GES plus ONIS) con- 
tinues outgoing through the DIG (Fig. [4) and reaches a 
large critical internuclear distance (Fig. 4] , Fig. [g] ,and 
Fig. [t] (d)) at which the ionization rate is dramatically 
enhanced p!5 l [2QH22] . Therefore, D2" is ionized consider- 
ably but mainly through the low-energy band of the KER 
spectrum (Fig.|3|. During 8 to 12.5 cycles, a considerable 
population of D2", before reaching to the larger critical 
internuclear distances, is ionized and also the laser field 
becomes weak (Fig. 1,4, and 5). Therefore, after this 



time interval, the ionization channel for the negligible 
residual population of D2" becomes nearly closed and the 
remaining D2" ions will dissociate at longer times. 
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